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FOREWORD 


For  many  years,  the  United  States  Air  Force  has  supported  experimental  inves¬ 
tigations  into  the  chemistry  and  physics  of  upper  atmospheric  energy  dynamics. 
In  January  1983,  to  further  this  research,  the  Air  Force  Geophysics  Laboratory 
(AFGL)  awarded  contract  F19628-83-C-0056  to  the  Space  Dynamics  Laboratories 
(SDL)  at  Utah  State  University  in  Logan,  Utah.  This  contract  obligated  SDL  to 
design,  build  and  operate  certain  field  instruments  whose  goal  was  to  measure 
the  infrared  spectral  and  spatial  characteristics  of  upper  atmospheric  energy 
transfer,  and  to  analyze  the  data  these  instruments  recorded.  As  one  step  in 
satisfying  this  obligation,  SDL  designed  and  built  a  rocket-borne,  field- 
widened  interferometer  (RBFWI)  and,  in  April  1983,  supported  its  launch  aboard 
a  Sergeant  rocket  (A30.276)  from  Poker  Flat  Research  Range,  Alaska,  into  a 
Class  II  aurora.  The  interferometer  recorded  significant  auroral  emission 
data  in  the  infrared  spectral  region  between  2.0  and  7,5  urn,  from  which  re¬ 
searchers  hoped  to  recover  vertical  profiles,  time  histories,  densities  and 
rotational/vibrational  temperatures  of  certain  atmospheric  constituents. 

Because  of  instrument  effects  inherent  in  the  RBFWI,  however,  ordinary  Fourier 
transform  analysis  of  RBFWI  interferograms  cannot  produce  spectra  of  suffi¬ 
cient  resolution  to  satisfactorily  recover  the  desired  information  with  stan¬ 
dard  spectral  inversion  processes.  To  overcome  this  problem,  SDL  sought  the 
services  of  Atmospheric  Radiation  Consul tants,  Inc.  (ARC),  of  Acton,  MA, 
awarding  them  subcontract  84-051.  Under  this  subcontract,  ARC  adapted  two 
computer  techniques  to  increase  the  resolution  of  the  RBFWI  spectra:  a  linear 
least-squares  analysis  and  a  nonlinear  deconvolution  algorithm.  This  publica¬ 
tion  describes  the  two  methods  and  summarizes  the  results  of  their  application 
to  the  RBFWI  data.  Unless  otherwise  noted,  all  the  work  reported  herein  was 
performed  by  the  authors  and  ARC,  Inc. 

The  authors  are  grateful  for  the  generous  support  and  encouragement  of  this 
project  by  Drs.  Ken  Schwartz  and  Peter  Lunn  of  the  Defense  Nuclear  Agency 
(DNA)  and  Drs,  Don  Ball  and  Ted  Cress  of  AFSOR. 
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SECTION  1 


INTRODUCTION  AND  SUMMARY 

The  Field-Widened  Interferometer  (FWI)  is  a  cryogenically-cooled , 
rapid-scan  Michelson  interferometer  spectrometer  with  special  field  com¬ 
pensation,  designed  to  measure  auroral  emissions  between  2  and  7.5  pm 
wavelength  at  ~1.1  cm  ^  resolution  (unapodized)  .  It  was  launched  aboard 
a  Sergeant  rocket  in  April,  1983  from  Poker  Flat,  Alaska  into  an  IBC 
Class  II  aurora.  In  this  experiment,  sponsored  by  the  Defense  Nuclear 
Agency  and  the  Air  Force  Geophysics  Laboratory,  the  uplooking  instrument 
obtained  interferograms  every  1.5  seconds  from  86  km  altitude  up  to 
139  km  (flight  apogee) ,  and  down  to  60  km  altitude  on  descent.  Details 
of  the  instrumentation  and  experiment,  and  examples  of  the  excellent  spec¬ 
tra  obtained,  are  given  by  Steed,  et  al  (1983). 

It  is  planned  to  analyze  spectra  from  the  FWI  experiment  by  an  inver¬ 
sion  process  that  will  recover  vertical  profiles  and  possibly  time  histo¬ 
ries  of  the  translational  temperature,  CO^  and  NO  excited  state  densities 
and  the  CO^  vibrational  temperature.  The  method  that  will  be  used,  a 
modification  of  the  limb  inversion  technique  developed  by  Zachor  and  Sharma 
(1985) ,  requires  the  use  of  individual  lines  (or  narrow  spectral  regions) 
of  differing  opacities,  including  optically  thin  lines.  The  FWI  spectra 
of  1,1  cm  ^  resolution  do  not  provide  adequate  spectral  isolation  because 
of  the  strong  sidelobes  in  the  apparatus  function  corresponding  to  this 
resolution.  Suppression  of  the  sidelobes  by  apodization  of  the  interfero¬ 
grams  degrades  the  resolution  to  '^2.2  cm  which  is  not  sufficient  to 
resolve  the  CO^  lines.  Thus,  the  inversion  technique  cannot  be  imple¬ 
mented  until  the  existing  FWI  spectra  are  a.)  analyzed  to  obtain  estimates 
of  individual  line  intensities  or  b.)  are  recomputed  by  a  method  that 
yields  higher  spectral  resolution  (i.e.,  higher  than  is  provided  by  Fourier 
analysis  of  the  interferograms) . 
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The  primary  purpose  of  the  present  study  was  to  set  in  place  proce¬ 
dures  and  computer  codes  for  performing  both  types  of  spectral  estimation. 
Straightforward  linear  least-squares  methods  were  adapted  for  the  deter¬ 
mination  of  line  intensities.  For  resolution  enhancement  we  selected  the 
constrained  nonlinear  deconvolution  algorithm  of  Jansson  (1984) .  A  second 
objective  was  to  determine,  by  these  techniques,  whether  ionized  nitric 
oxide  (NO^)  was  present  in  the  auroral  event  observed  by  the  FWI. 

Sections  2  and  3  describe  the  methods  and  results  of  their  application 
to  FWI  data  for  two  spectral  regions:  the  2300-2380  cm  ^  region  which 
spans  the  strong  band  of  CO^,  and  the  2220-2300  cm  ^  region  in  which  No"*^ 
would  probably  be  the  dominant  radiator,  if  present.  The  methods  are 
fairly  general  and  can  be  applied  to  data  that  will  be  obtained  in  forth¬ 
coming  AFGL/DNA  experiments,  including  HIRAM,  SPIRIT  and  CIRRIS  lA. 

The  figures  below  are  examples  of  results  obtained  for  the 
2300-2380  cm  ^  region  from  FWI  scan  10,  measured  at  ^'^101  km  altitude. 

The  first  figure  shows,  as  a  stick  spectrum,  the  intensities  of  the  CO^ 
lines  inferred  by  least-squares  with  the  constraint  that  the  intensities 
follow  a  smooth  band  contour.  This  display  is,  for  all  intents 
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SEGMENT 

(Fig.  18) 


and  purposes,  an  infinitely  resolved  spectrum  of  the  source  (even  the 
Doppler  line  shape  has  been  eliminated) .  The  least-squares  approach 
requires  knowledge  or  recognition  of  the  radiators  present  and  the  posi¬ 
tions  of  their  spectral  lines.  When  this  prior  information  is  available, 
the  method  is  both  the  simplest  and  most  effective  of  the  various  tech¬ 
niques  available  for  "superresolved"  spectral  recovery.  We  describe  in 
this  report  a  least-squares  method  that  is  applicable  even  when  overlapping 
bands  result  in  coincident  and  near-coincident  lines. 

The  second  figure  shows  the  spectrum  obtained  by  conventional  FFT 
analysis  of  the  triangularly  apodized  scan  10  interferogram,  and  a  section 
of  this  spectrum  after  deconvolution  by  the  Jansson  algorithm.  The 
apodization  is  necessary  because  the  Jansson  method  is  effective  only  if 
the  apparatus  function  is  positive  definite.  The  full  width  at  half  maxi¬ 
mum  (Ft'JHM)  of  the  lines  in  the  deconvolved  spectrum  is  ~0.5  cm  which 
is  4.4  times  better  than  the  '^2.2  cm  ^  resolution  of  the  starting  spectrum 
(or  2.2  times  better  than  the  1.1  cm  ^  resolution  for  the  unapodized  case). 
Moreover,  the  positions  of  these  CO^  lines  are  correct  to  within  ~0.2  cm  ^ 
and  there  are  no  sidelobes.  The  achieved  resolution  enhancement  is  much 
less  in  the  R  branch  of  the  CO2  band  where  the  lines  are  more  closely 
spaced.  Included  in  Section  3  is  a  spectrum  recovered  for  scan  10  by  the 
Maximum  Entropy  Method  (MEM) .  The  MEM  gave  somewhat  greater  resolution 
enhancement  but  did  not  accurately  reproduce  the  CO^  line  positions. 

Five  lines  of  the  NO^  1-0  band  (P4  through  P8)  and  four  lines  of  the 
2-1  band  (RO  through  R3)  fall  between  2310  cm  ^  and  2330  cm  but  none  of 
these  appear  in  the  deconvolved  spectrum  segment.  Similarly,  none  of  the 
many  lines  of  the  2-1  through  6-5  bands  that  occur  between  2220  and  2300  cm 
were  found  in  tt  j  deconvolved  FWI  spectra  for  this  region  (with  just  a  few 
exceptions  noted  in  Section  3.4) .  It  is  concluded  that  NO^  was  not  observed 
in  the  FWI  experiment.  There  are  many  emission  features  between  2220  and 
2300  cm  ^  in  the  deconvolved  spectra,  but  these  could  not  be  identified  with 
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any  common  atmospheric  species  such  as  C  O2  ,  N^O  and  that  radiate  in 


this  region.  The  emission  features  also  change  markedly  from  scan  to  scan, 
indicating  that  the  observed  source  was  changing  rapidly,  or  perhaps  that 
these  features  are  instrumental  artifacts. 

The  two  analysis  methods  described  and  implemented  in  this  study 
require  knowledge  of  the  instrumental  line  shape  or  apparatus  function. 

Since  it  was  not  determined  during  calibration,  we  assumed  it  to  be  the 
usual  sine  function  (for  the  unapodized  case).  It  is  recommended  that  the 
apparatus  function  be  measured  in  future  programs. 

The  FWI  experiment  included  ground-based  observations  of  the  auroral 

o  o 

event,  on-board  photometer  measurements  at  3914A  and  5577A,  and  on-board 
(in  situ)  measurements  of  oxygen  atom  density.  These  data  were  not  avail¬ 
able  to  us  and  hence  could  not  be  used  in  interpreting  the  measured  radiance 
spectra. 

In  summary,  the  reported  effort  demonstrates  the  remarkable  super¬ 
resolution  capaibilities  of  least  squares  analysis  and  t!.  '  Jansson  deconvolu¬ 
tion  method.  Application  of  these  methods  to  the  FWI  data  base  leaves  un¬ 
resolved  an  important  question,  apart  from  the  identification  of  a  radiator (s) 
in  the  2220-2230  cm  ^  region:  Why  do  the  CO^  lines  in  the  red  wing  of  the 
P-branch  seem  alternately  too  strong  or  too  weak,  whether  resolved  by  least- 
squares  (without  a  smoothness  constraint)  or  by  the  Jansson  method?  This 
behavior  is  seen  in  the  second  figure  atxjve.  Vibrationally  excited  CO^  is 
produced  by  energy  transfer  during  collisions  with  N2.  Since  the  radiative 
lifetime  of  the  mode  is  so  short  (~ 2  msec),  the  excited  CO^  molecules  at 
these  altitudes  (~  100  km)  are  probably  not  thermalized  before  radiating. 

The  observed  rotational  distribution  is  that  of  the  nascent  molecules  which 
may  not  be  Maxwell-Boltzmann,  The  observed  distribution  is  intriguing  but 
as  yet  we  have  no  explanation  for  it. 
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SECTION  2 


LINEAR  LEAST-SQUARES  ANALYSES  OF  FWI  SPECTRA 

2 . 1  Methods  of  Analysis 

Two  methods  will  be  described  that  estimate  the  intensities  of  lines 
in  a  measured  spectrum.  It  is  assumed  that  the  positions  of  lines  that  may 
be  present  are  known  to  high  accuracy,  and  that  the  apparatus  function  of 
the  spectrometer  is  known.  One  of  the  methods  solves  directly  for  the 
individual  intensities.  The  other  solves  for  a  set  of  coefficients  that 
define  the  band  spectral  contour,  or  several  contours  if  it  is  known  that 
there  are  overlapping  bands  in  the  spectrum.  These  coefficients  and  the 
line  positions  can  then  be  used  to  estimate  the  line  intensities. 

The  number  of  solutions  sought  in  either  method  is  always  smaller 
(usually  much  smaller)  than  the  number  of  spectrum  points  available.  That 
is,  the  solutions  give  best  estimates  in  the  overdetermi  led  least-squares 
sense.  The  methods  are  linear  and  are  therefore  easily  implemented  on  the 
computer. 


2.1.1  Estimation  of  Individual  Line  Intensities 


Let  L^j  denote  the  spectral  radiance  that  the  instrument  would  measure 
at  wavenumber  for  a  single,  very  narrow  line  source  of  wavenumber 
In  other  words,  L^^  is  a  shift-variant  apparatus  function  if  it  depends  on 
both  and  v^,  as  is  the  case  for  the  FWI,  In  the  more  usual  shift- 
invariant  case,  L^j  depends  on  only  Except  when  specifically  de¬ 

fined,  L^j  can  represent  either  case:  this  will  be  true  also  in  the 
equations  of  Section  3  that  define  the  Jansson  deconvolution  algorithm. 


The  source  emission  lines  in  the  FWI  experiment  are  much  narrower  than 
the  apparatus  function,  and  for  our  purposes  can  be  regarded  as  an  array 
of  delta  functions,  i.e.,  as  a  "stick  spectrum".  The  corresponding  measured 
spectrum  can  be  written 
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Given  the  measured  spectrum  (column  vector  S)  and  the  line  shape  matrix  L, 
one  can  obtain  an  estimate  a'  of  the  true  line  intensity  vector  a.  A  well- 
known  least-squares  solution  (one  that  minimizes  the  variance  between  S 
and  La')  is 

a'  =  +  alj  l'^S  ,  (3) 

T  -1 

where  (  )  and  (  )  denote  matrix  transpose  and  inverse,  and  I  is  the 

N  by  N  identity  matrix.  The  term  al  has  the  effect  of  suppressing  contribu- 

T 

tions  to  a  by  eigenvectors  of  L  L  with  sma) 1  eigenvalues,  i.e.,  of  suppres¬ 
sing  possible  instability  arising  from  noise  in  S.  An  appropriate  value 
for  a  if  the  system  noise-equivalent  spectral  radiance  (NESR)  is  known  is 

2 

Cl  =  (NESR)  /variance  of  S  ,  (4) 


Equation  3  is  useful  for  estimating  individual  line  intensities  when 
the  lines  are  well  separated.  Lines  may  be  coincident  or  nearly  so  when 
the  spectral  region  under  consideration  contains  two  or  more  overlapping 
bands  (bands  of  different  species  or  different  subbands  of  a  single  species) . 
Obviously,  lines  that  are  exactly  coincident,  must  be  treated  as  a  single 
line:  one  can  estimate  only  their  summed  intensities.  Equation  (3)  will 
give  poor  results  when  line  separations  occur  that  are  small  compared  to 
the  instrument  average  spectral  resolution  Av. 
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2.1.2  Estimation  of  Band  Contour  Parameters 

Band  "contour"  refers  here  to  the  smooth  function  or  envelope  describ¬ 
ing  the  variation  of  intensity  from  line  to  line  within  a  molecular  band. 
Spectral  lines  that  are  nearly  coincident  will  usually  belong  to  different 
bands,  and  it  is  unlikely  that  the  different  band  contours  will  be  coinci¬ 
dent  (or  proportional  to  one  another)  over  an  appreciable  fraction  of  the 
band  widths.  These  properties  of  overlapping  bands  and  knowledge  of  the 
line  positions  make  it  possible  to  solve  for  the  band  contours,  and  then, 
as  a  second  step,  for  the  individual  line  intensities.  If  the  contours  are 
represented  as  weighted  sums  of  basis  functions,  the  problem  remains  linear 
and  the  initial  solution  gives  estimates  of  the  weights.  This  solution  can 
be  cast  in  the  form  of  Eq.  (3) ,  as  will  be  demonstrated. 

For  later  convenience  we  change  the  summation  index  in  Eq.  (1)  from  j 
to  k  to  obtain 


S(v. ) 


^ik  ^ 


Again,  S(v.)  are  the  measured  spectrum  points.  The  a,  are  the  unknown 
intensities  of  lines  at  the  known  positions  If  the  lines  belong  to  a 

single  band,  the  a^^  can  be  regarded  as  samples  of  a  continuous  contour  or 
envelope  intensity  function  E(v),  which,  in  turn,  can  be  represented  by  a 
weighted  sum  of  continuous  basis  functions  Fj(v): 

\  =  E^^k’  -  S  ''j 

j 

where  the  c^  are  the  weights  (the  new  unknowns) .  The  result  of  substituting 
Eq.  (6)  into  Eq.  (5)  and  changing  the  order  of  summation  is 


S(v.) 


L . ,  F  .  (v,  )  C  . 
ik  D  k  3 


This  is  the  same  as 


-  E  hj 


Cj  f  i  \ f  ^  f 


.  =  >  L.,  F .  (v.  ) 

ij  /  j  ik  j  k 


Equation  (8)  has  the  same  form  as  Eq.  (1).  Thus,  a  least-squares  estimate 
of  the  N  weights  is  given  by  Eq.  (3)  after  appropriate  changes  in 
notation : 


The  matrix  element  ^ ,  in  physical  terms,  is  the  total  spectral 
radiance  at  due  to  all  source  lines  of  a  given  band,  but  whose  intensi¬ 
ties  include  only  the  jth  component  (basis  function)  of  the  band  contour. 
The  contribution  of  a  typical  line,  which  corresponds  to  a  single  term  in 
the  right  hand  side  of  Eq.  (9),  is  pictured  in  Figure  1. 

The  application  of  these  results  consists  of  the  following  steps: 
First,  choose  a  set  of  basis  functions  Fj(v)  capable  of  representing  a 
typical  band  contour.  Then,  from  the  given  line  positions  and  matrix 
elements  L^j^,  compute  matrix  G  using  Eq.  (9).  Solve  for  vector  c'  using 
Eq.  (10) ,  and  finally  for  the  line  intensities  using  Eq.  (6) . 
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The  discussion  of  the  method  thus  far  implies  that  the  lines  belong  to 
a  single  band,  i.e.,  that  there  are  no  additional,  overlapping  bands. 
Indeed,  the  present  method  may  be  preferable  to  the  one  described  in 
Section  2.1.1  in  the  non-overlapping  case  since  it  (the  present  method) 
imposes  the  constraint  that  the  estimated  line  intensities  follow  a  smooth 
contour.  The  choice  would  depend  on  available  signal-to-noise  and  the 
f lexib’ lity  of  the  selected  basis  functions;  a  smoothness  overconstraint 
is  obviously  undesirable. 


The  method  of  solution  described  in  this  section  is  made  applicable  to 

the  case  of  overlapping  bands  by  generalizing  the  definitions  of  matrices 

L  and  G.  Note  that  subscript  k  refers  always  to  a  particular  line  and 

thus  implies  a  particular  band.  For  example,  k  =  1  to  25  might  represent 

the  lines  of  band  A  and  k  =  2G  to  100  the  lines  of  band  B,  Similarly,  j 

refers  to  a  particular  basis  function  associated  with  a  particular  band; 

j  =  1  to  3  might  represent  the  basis  functions  for  band  A  and  j  =  4  a  .'^ijigle 

basis  function  used  for  band  B.  In  other  words,  multiple  bands  may  be 

accommodated  by  one  or  more  partitions  along  the  columns  of  L  and  G.  The 

column  partitions  in  G  correspond  to  row  partitions  in  the  column  vectors 

c  and  c' .  Figure  2  shows  the  partitions  in  G,  c  and  L  for  the  present 

example,  as  well  as  the  corresponding  column/row  partitions  in  F  =  [f  ] , 

representing  the  basis  function  samples  F.(v  ),  The  figure  also  illustrates 

3  ^ 

the  evaluation  of  Eqs.  (8)  and  (9) .  Note  that  certain  combinations  of  j 
and  k  are  not  allowed  (are  zero)  in  matrix  F,  Thus,  the  summation  indices 
in  Eqs.  (6)  and  (9) ,  which  both  involve  F,  are  restricted  to  certain 
values.  Equations  (8)  and  (10)  are  evaluated  in  the  normal  manner;  i.e., 
the  partitions  serve  merely  to  label  groups  of  vector  or  matrix  elements. 

The  labels  are  needed  in  Eq.  (10)  only  for  the  solution  vector  c'  so  that 
its  elements  can  be  assigned  to  the  proper  band,  and  so  that  the  estimated 
line  intensities  for  each  band  can  then  be  correctly  obtained  via  Eq,  (6) . 


We  conclude  this  subsection  with  a  general  description  of  the  basis 
functions  that  were  used  in  analyzing  FWI  spectra  Note  that  the  equations 
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require  only  the  values  F.(v  ),  equal  to  the  basis  functions  evaluated  at 
known  line  frequencies  A  particular  implies  a  particular  band,  a 

branch  of  that  band  (P,  Q  or  R) ,  and  a  rotational  quantum  number  J.  The 
contours  of  P  branches  were  represented  by  the  bases 


F.  {^)  )  =  J  exp[-a.  J(J+1)]  , 

1  ^  1 


where  J  has  the  value  implied  by  and  the  are  selected  from  the  set 
of  values 


J.  (2  J.  +  1) 
1  1 


(lla) 


determined  by  the  set  .  Equation  (11)  approximates  the  theoretical 
distribution  of  line  intensities  for  a  typical  P-branch  for  a  given  rota¬ 
tional  temiierature  (proportional  to  1/a  J  .  It  increases  initially  with  J 
(or  corresponding  distance  between  and  the  band  origin) ,  reaches  a 
maximum,  and  then  falls  off  with  a  long,  exponential-like  tail.  The  maxi¬ 
mum  occurs  at  J  =  J j .  One  might,  for  example,  choose  to  use  three  basis 
functions  whose  aVs  are  given  by  Eq.  (lla)  with  =  6 ,  10  and  14.  A 
weighted  sum  of  these  basis  functions,  representing  the  P  branch  contour, 
will  peak  anywhere  between  approximately  6  and  14  depending  on  the  relative 
weights . 

Similarly,  the  contours  of  R  branches  were  represented  by  the  bases 


F.(v  )  =  (J  +  1)  exp[-a.J(J  +  1)  ]  . 

1  -K  J 


The  a  ■  were  selected  from 
1 


+  1)  (2  J 


.tl)l-^ 

J 


(12a) 
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to  obtain  maxima  at  the  J  =  J . . 

3 

A  single  set  of  J^'s  was  used  for  the  P  and  R  branches  of  the  same 
band,  but  separate  solutions  (basis  function  weights)  were  obtained  for 
the  two  branches.  There  were  no  strong  Q  branches  in  the  FWI  spectra 
analyzed  by  this  method. 


2 . 2  Apparatus  Function  for  the  FWI 

The  shift-variant  apparatus  function  L. .  must  be  known  in  order  to 

ID 

apply  the  least  squares  methods  described  above  or  to  apply  the  deconvolu¬ 
tion  technique  described  in  Section  3.  Recall  that  L, .  =  L(v.,v.)  is 

ID  1  D 

defined  as  the  spectral  radiance  at  wavenumber  due  to  a  discrete  line 
source  of  wavenumber  u^.  The  interferograms  obtained  by  the  FWI  consist 
of  a  set  of  samples  corresponding  to  equal  increments  in  the  lateral  dis¬ 
placement  of  a  field-widening  wedge,  or,  equivalently,  to  equal  incre¬ 
ments  in  optical  retardation  x  at  some  reference  wavent  her.  Taking  the 
Fourier  transform  of  what(!ver  truncation/apodization  function  is  applied  to 
the  interferogram  yields  a  basic  shift-invariant  apparatus  function 
L(a^,aj),  where  o  is  the  variable  conjugate  to  the  sampling  retardation 
variable  x.  Conservation  of  energy  (line  intensities)  requires  that 


L(v.,v.)  dv .  -  L(o.,o.)  da.  , 

'  i'  3'  1  '  1'  1  ' 


i.e.,  the  integrals  over  the  two  sides  must  be  equal.  Note  that  v  is  some 
(nonlinear)  function  of  o  (or  o  a  function  of  v)  because  the  index  of 
refraction  of  the  w'^dge  m.iterial  varies  with  wavelength.  Thus,  if  L{o^,Oj) 
is  known,  the  desired  ap}' \ratus  function  =  L(v^,Vj)  can  be  obtained 


=  l[o(v^),  o(v^)]  (do/dv)^^^ 
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The  function  o (v)  and  its  derivative  do/dv  were  determined  by  USU  during 
calibration  of  the  FWI  and  have  been  provided  as  FORTRAN  function  subrou¬ 
tines  GFUNC(v)  and  GDERV(v),  respectively: 

a(v)  E  gfi;nc(v)-f 

da/dv  E  GDERV(v)-F  (15) 

F  E  5  X  632800  x  65536  x  10“^® 

-1 

V  given  in  cm 


The  two  subroutines  are  listed  below. 


FUNCTION  GFUNC(W) 

C 

C*  USED  IN  CALCULATING  THE  DISTORTED  FWI  LINESHAF'E 
C 

PARAMETER(  COEF 1  +  .  1 0252946  ) 

PARAMETER (  C0EF2=-. 15787537E5  ) 

PARAMETER <  C0EF3=-. 21755914  ElO  ) 

PARAMETER  <  C0EF4=  +  .  34318074E--1  1  ) 

WW  =  W*W 

GFUNC  =  W»<C0EFJ  +  <C0EF2  +  C0EF3/WW) /WW  +  C0EF4<tWW) 
RETURN 
END 
C 
C 
C 
C 
C 

FUNCTION  GDERV(W) 

PARAMETER <  COEF 1 =+. 1 0252946  ) 

PARAMETER<  C0EF2=-. 1 5707537E5  ) 

PARAMETER (  C0EF3=-. 2 1 7559 1 4E 1 O  ) 

PARAMETER <  C0EF4=+ . 343 1 8074E- 1 1  ) 

WW  =  W*W 

GDERV  =  COEFl  -  (C0EF2  +  3 .  *C0EF3/WW ) /WW  3.  *C0EF4)|!WW 

RETURN 

END 


The  least-squares  analysis  methods  were  applied  to  spectra  computed 
by  USU  from  unapodi/.ed  interferograms  of  length  64K  =  65536.  For  this 
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case  the  apparatus  function  in  the  c-space  is 


simi  [o  (v  . )  -o  (v  . )  ]  /Aa 

l[o  (v. )  ,0  (v  . )  ]  =  - - - - -  (16) 

(o  (v^)  -a  (v  j )  ] 

with  Aa  =  2.0.  The  correspondinq  wavenumber  resolution  Av  =  Ao‘(do/dv) 
where  the  derivative  is  given  by  Eq.  (15),  is  -^l.O  cm  ^  (~  0.91  cm  ^  to 
first  sine  zero  and  ~1.1  cm  ^  FWHM)  . 

The  FWI  spectra  in  the  region  2200-2320  cm  ^  change  marlcedly  from  scan 
to  scan  over  certain  portions  of  the  flight,  indicating  that  the  observed 
column  thickness  of  the  radiating  species  (NO+?)  is  changing  due  to  the 
motion  of  the  rocket,  and/or  that  local  properties  within  the  column  above 
the  rocket  are  changing  rapidly  with  time.  In  order  to  represent  the 
effects  of  a  "changing  source"  in  a  simplified  parametric  manner,  we  assume 
that  all  spectral  components  of  the  source  change  accord  ng  to  a  single 
multiplicative  function  of  time.  This  is  the  same  as  saying  the  effects 
can  be  represented  by  applying  a  single  apodization  function  to  the 
interferogram.  We  assume  that  this  function  is  an  exponential  and  charac¬ 
terize  it  by  its  value  c  at  maximum  retardation,  which  is  also  the  ratio 
of  apparent  source  intensity  at  maximum  retardation  to  that  at  zero 
retardation.  This  parameter  may  be  varied  to  obtain  maximum  agreement 
between  the  measured  spectrum  and  one  reconstructed  from  the  inferred  line 
intensities.  Note  that  c  can  be  greater  than  unity  for  source  intensity 
decreasing  with  time,  since  the  maximum  retardation  signal  is  measured 
before  the  zero  retardation  signal  during  even  numbered  scans. 

For  the  exponentially  varying  source,  the  basic  sine  function  instru¬ 
ment  line  shape  of  Eg.  (16)  is  replaced  by 

(1/Ao)  I  In  c  [c  cos  q  -  l]  +  c  q  sin  q  |  (17) 

[  (In  c)  ^  +  q^] 
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where  q  =  it  [a  (v^) -o  (v  j)  ]/Ao . 
for  c  1. 


Equation  (17)  reduces  to  the  sine  function 


2.3  Application  to  the  Band  (2300-2380  cm  ^  region) 

The  two  methods  described  above  were  applied  to  the  spectral  region 

2300-2380  cm  ^  of  FWI  scan  No.  10,  which  was  obtained  at  101  km  altitude 

during  the  upleg  portion  of  the  flight.  There  are  some  half-dozen  over- 
1 2  16 

lapping  bands  of  C  in  this  region,  and  it  was  necessary  to  decide 

which  CO2  bands/lines  to  include  in  the  analysis.  To  this  end,  theoretical 
or  "expected"  line  intensities  were  computed  using  AFGL's  code  NLTE  [see 
Sharma  et  al  (1983)],  which  had  just  been  modified  to  handle  the  zenith 
viewing  geometry  of  the  FWI  experiment  in  addition  to  the  limb  viewing 
geometry. 

12  16 

The  five  strongest  bands  of  C  for  this  region,  in  order  of 

decreasing  band  intensity  at  temperature  298K,  are  listed  in  Table  1.  The 
band  origins  and  lower  state  energies  are  taken  from  the  current  AFGL  line 
parameters  compilation  [Rothman,  et  al  (1983)].  The  NLTE  computations  are 
based  on  assumed  temperature,  pressure,  and  CO^  density  profiles,  and  an 
estimated  vibrational  temperature  profile  provided  by  J.  Kumer  of 
Lockheed  to  R.  Sharma  of  AFGL.  The  "observer"  altitude  was  taken  as  100  )cm 
The  stick  spectrum  of  Fig,  3  shows  the  computed  zenith-view  intensities  of 
all  lines  of  the  five  CO2  transitions  for  the  spectral  region  2300-2390  cm 
Lines  whose  apparent  intensities  are  less  than  the  vertical  resolution  of 
the  plot  (the  thickness  of  the  baseline  axis)  are  plotted  as  sticks  of 
height  equal  to  twice  this  resolution.  Only  the  two  strongest  bands  have 
definable  contours  in  the  figure. 


2.3.1  Individual  Line  Intensities 

It  was  decided  to  apply  the  method  of  Section  2.1.1  to  the  lines  of 
the  strongest  transition.  Note  that  the  maximum  line  intensity  for  the 


-18- 


Table  1.  The  Five  Strongest  C  O  Transitions  for  the 

2300-2380  cm  Region  [see  Rothman,  et  al  (1983)] 

Band  Origin  Band  Lower  State  Energy 

(cm  v'  yj"  (cm”^) 


2349.143 

00011 

00001 

0. 

2336.633 

01111 

01101 

667.379 

2324.141 

02211 

02201 

1335.129 

2327.433 

10012 

10002 

1285.4087 

2326.598 

10011 

10001 

1388.1847 

next  weaker  transition  (the  strongest  hot  band)  is  only  seven  percent  of 
that  for  the  included  lines.  Lines  and  data  outside  the  region  2300- 
2331.75  cm  ^  were  discarded;  this  resulted  in  a  total  of  N  =  52  lines  of 
CO2  (00011-00001)  and  M  =  328  measured  spectrum  values.  In  this  first 
analysis  we  also  ignored  the  possibility  that  overlapping  bands  of  NO^ 
were  present  in  the  measured  spectrum.  The  apodization  parameter  c  was 
set  to  unity,  corresponding  to  a  constant  source. 

The  theoretical  stick  spectrum  of  Fig.  3,  but  including  only  the  CO^ 
lines  whose  positions  were  used  in  the  analysis,  is  reproduced  in  the 
lower  part  of  Figure  4.  The  upper  part  of  the  figure  shows  the  USn  mea¬ 
sured  spectrum  for  scan  10.  Note  that  there  are  negative  amplitudes  due 
to  sidelobes  in  the  phase-corrected  USU  spectrum.  Note  also  that  the 
predicted  CO^  line  positions,  which  are  indicated  by  hash  marks,  do  not 
coincide  with  the  spectrum  peaks  except  near  v  0^2320  cm  ;  i.e.,  side- 

lobes  contribute  as  much  to  spectral  structure  as  the  central  lobes 
throughout  most  of  the  spectrum. 

Figure  5  shows  the  results  of  the  least-squares  analysis.  The  infer¬ 
red  line  intensities  are  displayed  as  a  stick  spectrum  in  the  lower  part  of 
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Figure  3.  C  0  line  intensities  computed  by  code  NLTE. 


ction:  Measured  spectrum  for  scan  10. 

ction:  Synthetic  (NLTE)  spectrum  of  CO.  00011-00001  lines. 


the  figure,  and  the  spectrum  synthesized  (reconstructed)  from  the  intensity 

★ 

estimates  is  shown  in  the  upper  part.  The  least-squares  process  has  mini¬ 
mized  the  variance  between  this  reconstructed  spectrum  and  the  measured  one. 
Figure  6  shows  the  inferred  line  intensities  again,  and  also  the  difference 
between  the  reconstructed  spectrum  and  the  measured  one.  The  rras  of  this 
difference  spectrum  is  only  four  percent  of  the  maximum  measured  spectral 
radiance;  the  maximum  difference,  which  occurs  near  v  =  2310  cm  ^ ,  is 
]3  percent  of  the  maximum  measured  spectral  radiance. 

The  maximum  spectral  radiance  in  the  measured  spectrum  (Fig.  4,  upper 
-11  2  -1 

curve)  is  3.9  x  10  W/cm  sr  cm  .  This  value  corresponds  to  a  maximum 

-13  2  -1 

signal-to-noise  ratio  of  445  based  on  the  NESR  of  2  x  10  W/cm  sr  cm 

(at  5  viin  wavelength)  quoted  by  Coo)c,  et  al  (1984)  for  the  FWI .  The  noise- 

to-signal  power,  denoted  by  a  in  Eqs.  (3)  and  (4),  was  found  to  be 
-5 

~ 5  X  10  ;  a  must  have  had  virtually  no  effect  in  the  solution  (Eq.  3) 

T 

since  the  diagonal  elements  of  L  L  were  or  order  10. 

The  inferred  line  intensities  (Fig.  5)  define  a  reasonably  smooth  band 
contour  at  wavenumbers  larger  than  ~  2325  cm  The  behavior  of  the  infer¬ 
red  intensities  below  this  wavenumber  might  at  first  be  thought  to  result 
from  near-coincidences  of  lines  in  the  00011-00001  and  01111-01101  bands. 

For  example,  Fig,  3  shows  that  the  P38,  36  and  34  lines  of  the  stronger 
band  occur  at  almost  the  same  wavenumbers  (~  2315,  2317  and  2319  cm  ^)  as 
three  lines  in  the  weaker  band.  However,  Fig.  5  shows  that  it  is  only  the 
P38  and  P34  lines  whose  inferred  intensities  are  greater  than  the  extrap¬ 
olated  band  contour.  This  very  same  behavior  is  observed  for  these  three 
lines  in  the  measured  sp>ectrum  after  its  resolution  has  been  enhanced  by 
spectral  deconvolution  (see  Section  3.2).  Moreover,  the  hot  band  lines 
are  conspicuously  absent  in  the  resolution-enhanced  spectrum. 

NOTE:  The  vertical  scales  for  the  stick  spectra  in  Figs.  4  and  5  are  not 

the'  same;  i.e.,  the  linos  have  been  nernalized  to  maximum  intensity. 
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For  the  reasons  indicated,  the  irregularities  in  the  inferred  00011- 
00001  intensities  were  attributed  to  overlapping  NO^  bands  rather  than 
the  strongest  overlapping  CO^  hot  band.  As  will  be  described,  we  attempted 
to  simultaneously  estimate  the  overlapping  contours  of  the  CO^  00011-00001 
band  and  several  NO^  bands  while  ignoring  possible  contributions  from  the 
CO.,  01111-01101  hot  band. 

Table  2  lists  the  CO^  line  positions,  their  inferred  intensities,  and 
the  intensities  computed  with  the  NLTE  code.  It  should  be  noted  that  the 
inferred  intensities  are  roughly  11  times  smaller  than  the  computed  ones. 
However,  the  contours  of  the  two  sets  are  very  similar  (compare  Figs.  4  and 
5,  lower  sections) . 

In  summary,  the  least-squares  method  for  estimating  individual  line 
intensities  has  yielded  excellent  results  for  the  band  of  CO^.  The  rms 
of  the  residual  error  is  only  four  percent,  which  is  rather  small  consider¬ 
ing  that  errors  would  be  expected  to  arise  from  imprecise  calibration  of  a 
nonlinear  wavenumber  scale  (i.e.,  from  fine-grained  spectral  variations  in 
index  of  refraction  of  the  field-widening/retarding  wedge) ,  from  urJenown 
effective  source  variations  during  the  scan,  and  from  instrument/detector 
noise.  Note  (Fig.  6)  that  if  the  region  below  ~  2330  cm  ^  and  the  region 
of  the  band  origin  are  excluded,  the  maximum  residual  error  is  roughly  five 
percent  and  the  rms  of  the  residual  is  considerably  smaller  than  the  four 
percent  quoted  above. 

2.3.2  Band  Contours 

12  16 

It  is  seen  from  Fig.  3  that  the  P  and  R  branches  of  the  C  00011- 

00001  band  have  pea)<s  near  J  —  12  or  14.  For  these  branches  we  used  the 
basis  functions  given  by  Kqs.  (11)  through  (12a),  with  =  2,  6,  10,  14, 

IS  and  22.  This  represents  a  total  of  six  basis  functions  for  each  branch 
or  12  for  the  entire  CO^  band.  As  a  test  to  determine  whether  this  was  an 
adequate  basis  set  we  performed  the  following  computations:  We  first 
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Table  2.  CO^  {00011-00001)  Line  Positions  {err.  )  and  Intensities  (W/cm  sr)  . 


convolved  the  stick  spectrum  representing  the  computed  (NLTE)  00011-00001 
line  intensities  with  the  shift-variant  FWI  apparatus  function  defined 

by  Eqs.  (14)  through  (16) .  This  theoretical  spectrum  was  then  analyzed  as 
data;  i.e.,  we  solved  for  the  vector  c'  of  weights  for  the  above-described 
basis  functions.  Using  the  solution  c'  we  then  reconstructed  the  spectrum, 
and  computed  the  rms  difference  between  it  and  the  original  theoretical 
spectrum.  This  rms  difference  was  only  0.2  percent  of  the  maximum  theoreti¬ 
cal  spectral  radiance.  Various  other  selections  for  the  J ^ ' s  resulted  in 
larger  rms  differences. 

The  spectral  region  2300-2380  cm  ^  also  includes  portions  of  the  fol¬ 
lowing  branches  of  NO^  bands: 

1- 0:  P  and  R 

2- 1:  P  and  R 

3- 2:  R 

4- 3:  R 

One  basis  function  was  used  for  each  of  these  six  branches.  Each  had  an  Qj 
corresponding  to  the  value  =  6,  based  on  computations  performed  by 

J.  Winick  of  AFGL.  The  analysis  included  the  contributions  of  lines 
located  within  "guardbands"  of  width  20  cm  ^  above  and  below  the  2300- 
2380  cm  ^  region  of  interest.  The  NO^  line  positions  required  in  the  anal¬ 
ysis  were  supplied  by  J.  Winnick. 

The  method  described  in  Section  2.1.2  was  used  to  estimate  the  scan  10 
CO^  band  contour  (two  branch  contours)  and  the  four  NO  band  contours  (six 
branch  contours  listed  above).  That  is,  we  solved  for  the  partitioned 
vector  c'  given  by  Eq.  (10) ,  which  determines  these  contours  and  the  corre¬ 
sponding  CO^  and  NO^  line  intensities.  Figure  7  shows  the  inferred  line 
intensities  plotted  as  stick  spectra.  The  upper  part  shows  only  the  CO^ 
intensities;  the  lower  part  shows  both  the  CO^  and  NO  intensities. 
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Table  3  lists  the  inferred  CO^  lino  intensities.  The  apodization  param¬ 
eter  c  is  unity  for  these  results. 

A  comparison  of  Tables  2  and  3  show  that  the  inferred  CO^  intensities 
have  not  changed  by  very  much,  but,  of  course,  the  new  estimates  follow  a 
very  smooth  contour.  Figure  7  shows  that  the  inferred  NO^  intensities  are 
very  low  in  comparison  to  those  for  CO^.  The  maximum  inferred  line  inten- 

4* 

sity  for  each  NO  band  relative  to  the  maximum  inferred  CO^  line  intensity 
is 

1- 0:  ~3.6  percent 

2- 1:  ^-<5.0  percent 

(~2.5  percent  for  R-branch  only) 

3- 2:  ~1.0  percent 

4- 3:  /n,6.0  percent  (occurs  near  2303  cm 

The  values  of  five  and  six  percent  for  the  2-1  and  4-3  L,..nds,  respectively, 
represent  lines  near  2300  cm  ^  where  the  residual  spectrum  error  is  high 
(~  10  percent  of  the  maximum  observed  spectral  radiance) .  The  inferred 
intensities  for  the  NO^  1-0  and  2-1  bands  are  plotted  on  normalized  scales 
in  Fig.  8. 

The  residual  spectrum  is  not  shown,  but  it  bears  a  strong  resemblance 
to  the  one  in  Fig.  6.  The  maximum  and  rms  residuals  relative  to  the  maxi¬ 
mum  observed  spectral  radiance  are  10  and  3.8  percent,  respectively,  com¬ 
pared  to  13  and  ~4  percent  for  the  method  represented  by  Fig.  6.  In  other 
words,  inclusion  of  NO^  has  not  effected  a  significant  reduction  in  the 
residual.  Based  on  tJieso  facts  and  the  very  small  ratio  of  inferred  NO 
intensities  to  CO^  intensities,  it  is  felt  that  the  present  results  are 
inconclusive  in  establishing  that  NO^  was  observed  by  tlie  FWI . 
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Figure  8.  Inferred  line  intensities  for  the  NO  1-0  and  2-1  bands. 

See  text  for  a  discussion  of  the  significance  of  these  results. 


2 . 4  Application  to  the  2220-2300  cm  ^  Region 

The  spectral  region  2220-2300  cm  ^  is  probably  optimum  for  the  detec¬ 
tion  of  NO^  in  the  FWI  spectra.  This  portion  of  the  spectrum  for  scan  10, 
unapodized,  was  analyzed  by  the  method  of  Section  2.1.2,  which  estimates 
overlapping  band  contours.  The  analysis  assumed  that  the  measured  emis¬ 
sions  are  due  to  the  following  NO^  bands/^ranches : 

1- 0:  P 

2- 1:  P 

3- 2:  P  and  R 

4- 3:  P  and  R 

5- 4:  R 

6- 5:  R 

Again,  one  basis  function  was  used  for  each  of  these  eight  NO^  branches. 

The  parameter  H  J^,  which  determines  the  location  of  the  basis  function 
peaks,  was  varied  between  6  and  10.  Only  lines  with  predicted  intensities 
greater  than  approximately  two  percent  of  the  maximum  predicted  line  inten¬ 
sity  in  each  branch  were  included;  J.  winick  provided  the  intensity  pre¬ 
dictions.  The  total  number  was  111  lines  of  NO^.  The  contributions  from 
CO^  lines  near  2300  cm  ^  or  higher  wavenumbers  were  assumed  to  be  negligible. 

The  measured  spectral  radiance  in  the  2220-2300  cm  ^  region  increases 
significantly  between  scans  9  and  10  and  between  scans  10  and  11;  i.e.,  it 

was  apparently  increasing  during  scan  10.  Since  the  optical  retardation 
decreased  with  time  during  the  scan,  it  is  appropriate  to  select  a  value 
less  than  unity  for  the  source-change  parameter  c.  This  parameter  was 
varied  from  0.005  through  unity  to  1.4. 

Figure  9  shows  the  USU  measured  spectrum  and  typical  spectra  recon¬ 
structed  from  the  least-squares  inferences.  The  agreement  between  the 
measured  and  reconstructed  spectra  is  obviously  very  poor.  For  c  =  0.3 
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SCAN  10,  MEASURED  SPECTRUM 


2228  2238  2248  2258  2268  2278  2288  2298  2308 

WAVENUMBER,  cm  ^ 


Figure  9.  Measured  spectrum  for  scan  10  and  two  spectra  reconstructed 
from  inferred  NO'*’  band  contours  and  line  positions. 


-33- 


the  rms  difference  is  approximately  22  percent  and  the  maximum  difference 

is  approximately  67  percent;  for  c  =  1.0  these  differences  are  27  percent 

and  86  percent.  The  smallest  rms  difference  obtained  for  the  entire  range 

of  variation  of  c  and  was  approximately  21  percent  Similarly,  the 

maximum  difference  could  not  be  reduced  below  approximately  57  percent. 

+ 

These  results  indicate  that  the  measured  spectrum  is  not  due  to  NO  . 


SECTION  3 


DECONVOLUTION  ANALYSTS  OE  FWl  SPECTRA 

3 . 1  Backgrou nd 

The  least-squares  methcxls  described  above  rely  on  knowledge  of  the 
particular  radiators  present  in  the  source  and  on  the  positions  of  their 
spectral  lines.  If  this  information  is  lacking,  or  if  errors  in  the  cali¬ 
brated  wavelength  scale  are  suspected,  these  methods  arc  not  applicable  or 
marginally  useful.  The  present  section  describes  and  applies  a  technique 
for  recovering  " superresolved "  spectra  in  Fourier  transform  spectroscopy , 
i.e.,  of  obtaining  spectral  resolution  higher  thaii  that  provided  by  conven¬ 
tional  EFT  analysis  of  inter ferog rams .  If  the  achieved  resolution  is  high 
enough,  the  spectra  nay  reveal  sequences  of  lines  (or  band  contours)  that 
are  recognized  as  belonging  to  a  particular  species;  resolved  (isolated) 
lines  can  be  integrated  to  obtain  their  intensities.  Wi*-h  these  possibili¬ 
ties  in  mind  we  coded  the  technique  and  applied  it  to  the  2220-2300  cm  ^ 
and  2300-2380  cm  ^  regions  of  several  FWI  spectra.  In  particular,  we 
wanted  to  find  out  which  species,  if  not  NO^,  is  (are)  responsible  for  the 
observed  emissions  between  2220  and  2300  cm 

Modern  methods  of  signal  recovery  and  spectral  estimation  include  the 
Maximum  Entropy  Method  (MEM) ,  the  related  Linear  Prediction  (LP)  and  Auto- 
Regressive  (AR)  techniques.  Maximum  Likelihood  (ML) ,  Minimum  Cross  Entropy 
(MCE) ,  and  constrained  nonlinear  deconvolution.  They  are  often  used  with 
the  objective  of  removing,  to  the  extent  possible,  the  blurring  effect  of 
the  measuring  instrument,  i.e.,  to  recover  fine  structure  in  a  measured  or 
inferred  spectrum  or  spatial  detail  in  an  image.  For  this  type  of  applica¬ 
tion  they  may  all  be  characterized  as  constrained  (linear  or  nonlinear) 
methods  of  deconvolution,  even  though  some  operate  on  transforms  of  the 
desired  signal.  The  methods  referred  to  have  beer,  demonstrated  to  achieve 
resolutions  substantially  greater  than  the  limits  implicit  in  the 
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"Blackman-Tukey "  approach  based  on  Fourier  analysis.  These  limits  are 
defined  by,  for  example,  the  width  of  an  entrance  slit,  length  of  an  inter- 
ferogram,  or  the  size  of  a  collecting  aperture  or  detector.  A  comprehensive 
summary  of  the  methods  available  for  deconvolution  is  given  by  Jansson 
(19R4)  . 

We  first  atti'mpted  to  recover  superresolved  FWI  spectra  in  the  region 
of  the  CO^  band  (2300-2380  cm  by  the  MEM.  The  result,  shown  in 
Fig.  10,  exhibits  supperresolution,  but  the  CO^  line  positions  are  not 
accurately  reproduced.  Specifically,  the  line  spacing  in  the  R-branch  is 
just  slightly  less  than  that  in  the  P-branch;  the  ratio  at  the  branch 
maxima  ought  to  be  ~  1.35/1.8  =  0.75.  Also,  the  gap  at  the  band  origin 
is  too  large.  The  wavenumber  scale  in  Fig.  10  is  not  calibrated.  However, 
calibration  would  not  produce  a  significant  change  in  the  ratio  of  line 
spacings  for  the  two  branches. 

There  are  many  examples  in  the  literature  of  the  successful  use  of  the 
Mi'M  to  obtain  superresolution  (with  correctly  reproduced  line  positions)  in 
the  case  of  spectra  with  ~ 10  or  fewer  lines  that  are  not  uniformly  spaced. 
There  have  not  been,  at  least  to  our  knwoledge,  any  similarly  successful 
demonstrations  for  spectra  containing  a  large  number  of  lines  (say,  80  or 
more)  with  pseudo-regular  spacing  like  the  CO^  band. 

After  this  experience  with  the  MEM  we  selected  the  Jansson  method  of 
constrained  nonlinear  deconvolution  for  use  in  analyzing  FWI  spectra; 
see  Jansson  (1984) .  The  latter  has  been  proven  to  work  well  for  arbitrarily 
complicated  spectra,  although  its  superresolution  capability  is  not  nearly 
as  dramatic  as  that  of  the  MEM  applied  to  simple  spectra. 

3 . 2  Jansson 's  Method 

Jansson 's  method  of  constrained  nonlinear  'leconvolution  is  an  out¬ 
growth  of  the  linear  ileconvolution  method  of  Van  Cittert.  The  nonlinearity 
of  the  Jansson  algorithm  is  the  result  of  a  constraint  that  suppresses 
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unphysical  neqative  sideloads  that  tend  to  be  generated  by  the  Van  Cittert 
algorithm  (even  when  the  apparatus  function  is  positive  definite) ;  see 
Jansson  (1984),  pp.  71-73.  This  constraint  results  in  a  greatly  improved 
capability  to  enhance  spectral  resolution. 

The  following  notations  will  be  used  in  describing  the  Jansson 
algorithm.: 


O. 


1 


I  . 
1 


L  .  .  (or  I, ,  . ) 

1]  1-] 


the  sought-after  "object"  spectrum,  i.e.,  the  true  source 
spectrum , 

the  "image"  or  observed  (degraded)  spectrum,  formerly 
denoted  S  (v^) , 

the  shift-variant  (or  invariant)  apparatus  function, 
superscript  denoting  the  iteration  cycle. 


The  algorithm  is 


+  R(of^)  {l.  -(L*o' 


(18) 


where  *  denotes  convolution,  and  R  is  a  "relaxation  function"  to  be 
defined  subsequently.  The  superscripts  on  the  involved  in  the  con¬ 
volution  are  not  defined,  to  indicate  that  the  "point-successive"  method 
of  evaluation  is  used.  That  is,  values  of  that  have  already  been  cal¬ 
culated  in  the  current  (k  + 1st)  cycle  are  used;  the  remaining  values  of 
represent  the  previous  (kth)  cycle.  Thus,  the  convolution  term  is 


(l  *  O 


(  ) 


E 


L., 

^3  3 


+ 


(19) 


r-; 


i 


This  expression  belies  the  simplicity  of  the  algorithm.  Note  that  the  O. 

(k)  ^ 

can  be  stored  in  a  single  array;  any  element  0^  is  recomputed  as 
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o|k+l)  ^  convolution  in  Fq.  (18)  uses  the  current  array, 

i.e.,  the  array  elements  corresponding  to  the  current  values  of  i  and  k. 
This  method  of  evaluation  has  been  found  to  be  more  rapidly  convergent 
than  the  "point-simultaneous"  method,  which  uses  the  convolution  of  L  and 
,  and  hence  requires  the  storage  of  both  and 

Equation  (18)  is  the  Van  Cittert  algorithm  if  the  relaxation  function 
R  is  set  to  unity.  In  this  case  the  difference  between  the  current  and 
previous  estimates  of  the  object  spectrum  is  equal  to  the  difference  be¬ 
tween  the  image  spectrum  I  and  the  current  "reconstructed  image"  spectrum 
L  *0^  ^ .  This  difference  is  zero  if  ^  =  I  or  when  the  process  has 
converged. 

Generally,  the  initial  guess  is  taken  as  I.  Suppose  the  true 

spectrum  (object)  is  a  single  isolated  line.  Then  will  be  broader 

than  I,  and  the  bracketed  term  in  Eq.  (18)  will  be  positive  near  the  line 

center  and  negative  in  the  line  wings  (assuming  the  app..  atus  function  L 

is  normalized,  positive  definite  and  monotonically  decreasing  away  from  the 

line  center,  like  a  Gaussian) .  Consequently,  0^^^  ,  evaluated  by  the  Van 

Cittert  algorithm,  will  be  bigger  than  0^*^^at  the  line  center,  and  will 

have  negative  "sidelobes".  To  suppress  this  undesirable  effect,  Jansson 

(k)  (k) 

chose  a  relaxation  function  R(0^  )  that  is  negative  when  is  negative. 

In  fact,  Jansson  used  the  relaxation  function  pictured  in  the  upper  part 

(k) 

of  Fig.  11,  which  is  negative  when  is  negative  or  greater  than  unity. 

Jansson  was  interested  in  deconvolving  absorption  spectra;  i.e.,  values 

(k) 

outside  the  range  0  0^  ^1  are  unphysical.  This  relaxation  function  is 

given  by 


"■0 

-  l/2\ 

). 

(20) 

/ 

L.  . 

X  1 

\ 

*  X  ’ 

i 

1 

where  is  an  adjustable  parameter  tliat  controls  the  "strength"  of  the 

relaxation  function. 


laxation  function  used  in  the  present  study. 


since  the  present  study  is  concerned  with  emission  spectra,  we  used 
the  inonotonically  increasing  relaxation  formula 

2  of  Vmax[o{''"^’]  (21) 

ii 

which  is  pictured  in  the  lower  part  of  Fig.  11.  Equations  (20)  and  (21) 

(k) 

give  approximately  the  same  value  when  is  less  than  one-half  its 

maximum  value.  In  our  application,  O  and  I  have  units  of  spectral 
radiance . 

The  Jansson  method  with  the  minimum  negativity  constraint  imposed 
by  the  relaxation  function  (20)  or  (21) ,  presupposes  a  positive  definite 
apparatus  function.  This  can  be  achieved  in  Fourier  transform  spectros¬ 
copy  by  appropriate  anodization  of  the  inter f erogram .  T)'ie  apodization 
will  result  in  an  immediate  loss  in  spectral  resolution  ^equivalent  to  a 
factor  of  two  or  more),  but  will  also  reduce  the  amplitude  of,  or  even 
eliminate,  the  sidelobes  of  the  apparatus  function.  We  elected  to  use 
triangular  apodization.  It  will  be  demonstrated  that  the  Jansson  algorithm 
can  recover  much  more  than  the  factor-of-two  initial  loss  in  resolution  and 
virtually  eliminate  the  residual  "artifacts",  i.e.,  the  sidelobes  of  the 
new  sine-squared  apparatus  function. 

We  did  not  actually  apply  triangular  apodization  to  the  FWI  interfero- 
grams,  but  rather  convolved  the  FWI  spectra  with  the  corresponding  shift- 
variant  apparatus  function.  This  apparatus  function  is  obtained  by  substi- 
tuting  [(sin  q/2)/(q/2)]  /(2Aa)  for  L(c^,Oj)  in  Eq.  (14);  i.e.,  this  sinc- 

sq'iared  function  replaces  the  sine  function  (sin  q/q)/Ao  given  by  Eq.  (16). 
By  performing  spectral  convolution  rather  than  interferogram  apodization, 
we  avoided  having  to  repeat  the  calibrations  (both  wavelength  scale  and 
spectral  responsivity)  and  interpolation  performed  by  USD. 


Recall  that  in  the  least-squares  analysis  of  FWI  specti a  we  allowed 
for  an  exponential  source  variation  characterized  by  the  parameter  c 
(Section  2.2).  To  represent  the  same  effect  in  the  present  deconvolution 
analysis,  the  apparatus  function  L  must  be  modified  to  include  the  source 
exponential  self-apodization:  L  is  given  by  Eq.  (14)  with 


L(a^,aj) 


e  inct^  cosq 


,  2c(ln  c)q  ,,  N 

+  - ^ - —  sin  q  -  (In  c)  ( 1+c)  - 


where 

N  H  (In  c)  ^  -  q^ 
D  =  (lnc)^  +  q^ 


(23) 


and  where  q  =  n (a (v^) -d (v ^ ) ] /Aa ,  as  before.  This  apparatus  function 

L(a.,a.)  =  L(q)  is  the  Fourier  transform  of  an  exponential  exp(|-ax|) 

^  ^  2 
times  a  triangle.  It  reduces  to  the  sine  (q/2)  function  described  above 

who  n  c  =  1 . 

Note  that  two  different  apparatus  functions  are  involved  in  the 
deconvolution  analysis.  The  first  is  convolved  with  the  USU-computed 
spectrum  to  produce  the  effect  of  triangular  apodization.  Any  effects  of 
a  varying  source  are  already  present  in  this  image  spectrum.  The  second 
apparatus  function,  defined  by  Fqs.  (14),  (22)  and  (23),  is  used  in  the 

deconvolution  algorithm,  Eq.  (18) .  It  accounts  for  the  triangular  apodiza¬ 
tion  and  the  possibility  of  a  varying  source  (c  ^  1). 

Jansson  points  out  the  importance  of  smoothing  the  original  data  and 
also  of  smoothing  the  result  of  each  iteration  cycle  to  prevent  the 
buildup  of  high-frequency  noiae;  see  Jansson  (1984)  and  wor)cs  cited  on 
pg.  106  of  this  reference.  I’liv.i,  ct  al  (1980)  emphasize  the  smoothing 
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requirement  as  well  as  the  need  for  high  signal-to-noise  in  the  original 
data  spectrum.  We  used  essentially  the  same  si-ioothing  procedures  as 
Jansson,  Like  Jansson,  we  used  a  smoothed  version  of  the  data  spectrum  as 
the  initial  guess  for  the  object  spectrum. 

The  FWI  spectra  provided  by  USU  consist  of  points  spaced  0.25  cm  ^ 

_  1 

apart.  As  mentioned  above,  the  spectral  resolution  is  approximately  1  cm  * 
unapodized,  or  ~2  cm  ^  apodized.  Tiie  0.25  cm  ^  spacing  is  not  adequate 
if  convolutions  are  performed  by  the  simple  rectangular  integration  implied 
by  Eq.  (19) .  We  therefore  interpolated  the  data  to  obtain  eight  times  as 
many  points.  The  interpolation  was  performed  as  part  of  the  initial  con¬ 
volution  that  effected  triangular  apodization;  i.e.,  the  convolution  lags 
were  multiples  of  0.25/8  cm 

The  computer  code  written  to  perform  the  Jansson  deconvolution  was 
applied,  during  the  final  debugging  phase/  to  the  case  in  which  the  image 
spectrum  is  due  to  a  source  consisting  of  just  a  few  del*- a-f unction  emission 
lines.  The  artificial  image  spectrum,  after  the  initial  convolution/ 
interpolation  step,  consisted  of  sine-squared  lines  whose  widths  correspond 
to  spectral  resolution  Av  cs- 2  cm  \  One  such  spectrum  is  shown  in  the  top 
panel  of  Fig.  12,  It  has  two  lines  separated  by  2.5  cm  ^  and  a  third, 
isolated  line.  The  middle  and  bottom  panels  show  the  deconvolved  spectrum 
after  ten  and  29  iterations,  respectively.  For  all  deconvolut,.ons  per¬ 
formed  in  this  study  the  lelaxation  strength  parameter  was  fixed  at  0.1 
for  the  first  half  of  the  itercations  and  at  0.2  for  the  second  i'.aif. 

Note  that  the  resolution  (based  on  FWHM)  is  enhanced  by  roughly  a 
factor  of  three  in  this  n'ise-free  case,  and  that  the  sidelobes  have  all 
but  disappeared  after  20  iterations.  Most  of  the  increased  resolution  is 
gained  in  the  first  10  cycles,  hut  the  relative  sidelobe  amplitude  continues 
to  'Increase  between  the  loth  and  the  20th  cycle.  It  should  be  noted  that 
the  maximum  in  each  sneetTum  is  plotted  to  the  same  height;  a  plot  in 
absolute  units  would  demonstrate  that  the  line  intensities  are  preserved. 
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Figure  13  shows  the  initial  spectrum  and  deconvolved  result  after  20 
iterations  when  the  two  closely  spaced  lines  are  1.25  cm  ^  apart,  i.e., 
are  unresolved  in  the  image  spectrian.  The  algorithm  succeeds  in  resolving 
the  two  lines,  but  does  not  eliminate  their  overlap  entirely,  despite  the 
fact  that  the  final  base  width  of  the  isolated  line  is  approximately 
1.25  cm 

3.3  Application  to  the  CO^  Band 

The  Jcinsson  method  was  used  to  deconvolve  the  spectral  region 
2290-2390  cm  ^  of  FWI  scans  10  and  6.  The  objectives  of  this  effort  were 

a. )  to  further  validate  the  method  and  computer  code,  based  on  the 

known  positio.ns  of  the  CO^  lines. 

b. )  to  test  the  scan-to-scan  reproducibility  of  the  FWI  spectra 

with  respect  to  line  positions,  and 

c.  )  to  determine  if  NO^  lines  of  the  1-0  and  2-1  bands  are  present 

between  the  CO2  lines 

The  spectrum  computed  by  USU  for  scan  10  from  a  64K  unapodized  inter- 
ferogram  was  already  shown  in  Fig.  4.  Convolution  of  this  spectrum  with 
the  sine-squared  function  that  effects  triangular  apodization  yields  the 
spectrum  shown  in  Fig.  14.  This  is  the  starting  point  (image  spectrum) 
for  the  deconvolution  algorithm.  It  represents  very  nearly  the  same  spec¬ 
tral  resolution  that  results  from  Hamming  apodization  of  a  64K 
interferogram.  Note  that  very  little  spectral  structure  remains  in  the 
R-branch  of  the  CO.,  band,  where  the  average  line  spacing  is  approximately 
1.35  cm. 

Figure  15  shows  the  deconvolved  spectrum  after  60  iterations.  For 
this  result  the  source  apodization  parameter  is  1.05,  representing  an 
assumed  five-percent  decrease  in  the  CO.,  column  thickness  during  scan  10 
(a  "reverse"  scan).  Various  other  values  up  to  1.4  were  tried,  but  these 
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deconvolving  the  spectrum  of  Fig.  14;  60  iterations 


resulted  in  less  resolution  enhancement.  Note  that  the  resolution  enhance¬ 
ment  is  very  good  in  the  red  wing  of  the  P-branch  where  the  line  spacing 
is  approximately  2  cm  ^ .  The  enhancement  is  relatively  poor  at  wavenumbers 
larger  than  ~2335  cm 

The  same  deconvolved  scan  10  spectrum  is  shown  on  an  expanded  wave- 
n’jmiier  scale  in  Fig.  16,  which  consists  of  five  panels  covering  the  spec¬ 
tral  region  2290-2390  cm  Figure  17  shows  the  corresponding  deconvolved 
spectrum  for  scan  6.  Note  that  the  plots  are  normalized;  the  maximum 
spectral  radiance  for  scan  6  is  roughly  twice  that  of  scan  10.  A  notice¬ 
able  feature  of  the  result  for  scan  6  is  the  saturated  spectral  radiances 
for  many  of  the  P  branch  lines,  indicating  that  the  source  is  optically 
thick  in  these  lines.  The  FWI  data  for  scan  6  was  obtained  at  ~ 95  km 
altitude,  whereas  scan  10  data  represents  ~ 101  km  altitude. 

Figure  18  reproduces  the  2310-2330  cm  ^  region  of  the  deconvolved 
spectrum  for  scan  10.  Several  features  are  noteworthy:  The  positions  of 
the  CO.,  lines,  especially  P42  through  P32,  are  in  excellent  agreement  with 
the  positions  given  in  the  AFGL  line  parameters  compilation,  which  are 
indicated  by  tl'.e  hash  marks.  Secondly,  there  are  no  lines  other  than  the 

expected  00011-00001)  lines.  In  particular,  lines  of  the  1-0  and 

+  * 

2-1  bands  of  NO  ,  which  overlap  the  CO^  band,  are  absent.  Note  also  that 

the  CO^  P42  through  P28  lines  seem  alternately  too  strong  or  too  weak.  The 
same  alternating  intensity  distribution  was  observed  in  the  results  of  the 
least-squares  analysis  of  scan  10  (see  Section  2.3.1  and  Fig.  5). 

Therefore  it  is  not  an  artificial  effect  of  the  Jansson  method.  Finally, 
the  FWUM  of  the  deconvolved  lines  is  ~0.4  to  0.5  cm  ^  compared  to 
~2.2  cm  ^  in  th.e  image  spectrum.  This  represents  a  resolution  enhancement 
factor  of  4.4  (or  ~2.2  if  the  final  resolution  is  compared  to  the  basic 
1.1  cm  ^  capability  of  the  FWI). 


*  + 

These  NO 


lines  have  the  positions  showmi 


in  Fig. 


8. 
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Figure  17.  Sair.f?  as  Fig.  16,  except  for 


scan  6. 
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The  2310-2330  cm  sccinents  of  fhe  doconvolveci  spectra  for  scans  10 
and  f>  are  conpared  in  Fiy.  10.  It  would  appear  that  the  wavenumber  scale 
of  scan  6  is  slighly  shifted  and  stretched  relative  to  the  scale  for  scan 
10.  Also,  the  resolution  enhancement  for  scan  6  is  noticeably  less  after 
the  first  few  lines  (P42  through  P38) .  This  difference  in  achieved  resolu¬ 
tion  is  possibly  related  to  the  existence  of  both  thin  and  thick  lines  in 

the  spectrum:  the  part  of  the  apparatus  function  due  to  source  variations 
may  be  quite  different  for  the  thin  and  thick  regions.  We  used  c  =  1.05 

and  60  iteration  cycles  in  deconvolving  scan  6;  c  was  not  varied.  Perhaps 

this  value  of  c  was  not  appropriate  for  scan  6  (or  for  the  self-absorbed 
lines  of  scan  6) . 

Unexplained  shifts  are  also  present  in  spectra  computed  by  USU  from 
64k,  Hamming-apodized  inter ferograms.  Four  such  spectra  representing  four 
consecutive  scans  are  shown  in  Fig.  20.  Spectrum  11  appears  to  be  shifted 
relative  to  spectrum  10. 


3 . 4  Application  to  the  2220-2300  cm  Region 


FWI  data  for  this  s[  ectral  region  were  deconvolved  in  the  same  manner 
as  for  the  2300-2380  cm  ^  region,  except  that  only  20  iterations  were  used 
(since  60  iterations  required  considerable  computer  time) ;  also  c  was 
fixed  at  0.95.  The  results  were  examined  for  NO^  lines  even  though  the 
least-squares  analysis  of  this  region  (Section  2.4)  and  the  deconvolutions 
described  above  (Section  3.3)  indicate  that  was  not  observed  by  the 


FWT . 


We  elected  to  use  scans  8,  6,  2  and  the  sum  (coaddition)  of  scans  2, 

3  and  4.  Scans  8  and  6  were  selected  because  they  have  little  or  none  of 
what  resembles  a  "continuum"  component  in  most  other  scans.  As  pointed  out 
by  Pliva,  et  al  (1980),  a  continuum  can  result  in  spurious,  weak  lines 
(similar  to  sidelobes)  adjacent  to  the  deconvolved  emission  lines. 

Scans  2,  3  and  4  were  used  because  they  were  tlie  basis  for  identification 
of  NO"^  by  USU. 
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Figure  20.  Spectra  computed  (by  USU)  from  Hamming-apodized  interferogram.  Note  the 
aooarent  wavenumber  scale  shift  between  scans  10  and  11. 


Figures  21  through  24  show  the  deconvolution  results  for  scans  8,  6, 

2  and  2+3+4,  respectively.  An  obvious  visual  feature  of  these  spectra 
is  their  overall  poor  correlation  with  one  another.  Indeed,  the  initial 
(image)  spectra  exhibit  a  similar  lack  of  correlation,  as  do  the  spectra 
computed  by  USU  with  Hamming  ai)odization .  There  are,  however,  a  few  lines 
that  appear  consistently  in  the  four  spectra:  the  ones  near  2233  cm 
2252  cm  2258  cm  ^  and  the  one  just  short  of  2290  cm 

Figure  25  shows  the  calculated  locations  of  NO^  lines  with  rotational 
quantum  numbers  up  to  ~  20  in  the  2220-2300  cm  ^  region.  Each  of  the  bands 
1-0  through  6-5  are  represented  by  a  displaced  set  of  hash  marks,  and  again 
by  shorter  hashes  just  above  the  wavenumber  axis.  J.  Winick  of  AFGL  pro¬ 
vided  the  line  positions  as  well  as  relative  line  strength  estimates  (not 
indicated  in  the  figure)  for  an  assumed  rotational  temperature.  For  this 
temperature,  the  strongest  lines  occur  near  J  -■  7  in  the  P  branches  and 
near  J  =  5  in  the  R  branches.  The  strongest  lines  of  each  band  are  indi¬ 
cated  by  the  emphasized  hash  marks. 

A  transparency  of  Fig.  25  was  overlaid  on  the  four  deconvolution 
results  (Figs.  21-24)  to  determine  if  the  latter  contain  a  significant 
number  of  the  stronger  NO^  lines.  The  predicted  positions  of  the  P5 
through  P9  lines  of  the  2-1  band  coincide  with  relative  maxima  in  the 
deconvolved  coadded  spectrum  (Fig.  24).  However,  there  are  few,  if  any, 
similar  correlations  in  other  parts  of  this  spectrum,  or  in  the  other  three 
cases,  even  if  one  allows  for  slight  spectral  shifts  or  scale  stretching. 
Our  conclusion  is  that  NO^  is  not  present  or  is  masked  by  a  stronger 
emitter(s) . 

other  common  atmospheric  species  that  radiate  in  the  2220-2300  cm  ^ 
region  are  N^O  and  ozone.  A  cursory  review  of  laboratory  data, 

solar  spectra  and  line  atlases  indicated  that  these  molecules  are  not 
responsible  for  the  observed  emissions. 
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Fi'^ure  21.  Deconvolved  spectrum  for  scan 


wavenumber,  cm 


WAVE^^JMBER,  cm 


WAVENUMBER,  cr. 


We  have  one  reservation  that  should  bo  stated  regarding  the  cioconvolu- 
tion  results  for  the  2220-2300  cm  ^  region.  The  signal-to-noise  for  this 
region  is  estimated  to  be  only  ^<45,  Pliva,  et  al  (19P0) ,  in  deconvolving 
Doppler-limited  spectra,  state  that  when  their  signal-to-noise  was  reduced 
to  100  they  could  "reliably  deconvolve  only  to  a  resolution  enhancement 
factor  of  2-2.5  before  objectionable  noise  begins  to  appear  on  the  de¬ 
convolved  curve".  Our  results  may  contain  artifacts  due  to  the  relatively 
low  signal-to-noise.  On  the  other  hand,  the  scan-to-scan  nonreproducibility 
of  the  image  spectra  suggests  that  instrumental  artifacts  may  dominate  the 
spectrum  in  low  signal-to-noise  regions. 
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